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ABSTRACT 

We describe a set of simulations of super-critical accretion onto a non-rotating supermassive 
BH. The accretion flow is radiation pressure dominated and takes the form of a geometri¬ 
cally thick disk with twin low-density funnels around the rotation axis. For accretion rates 
> lOMEdd. there is sufficient gas in the funnel to make this region optically thick. Radiation 
from the disk first fiows into the funnel, after which it accelerates the optically thick funnel 
gas along the axis. The resulting jet is baryon-loaded and has a terminal density-weighted 
velocity 0.3c. Much of the radiative luminosity is converted into kinetic energy by the time 
the escaping gas becomes optically thin. These jets are not powered by black hole rotation or 
magnetic driving, but purely by radiation. Their characteristic beaming angle is ~ 0.2 radians. 
For an observer viewing down the axis, the isotropic equivalent luminosity of total energy 
is as much as 10^^ ergs“^ for a 10^M© BH accreting at 10^ Eddington. Therefore, energeti¬ 
cally, the simulated jets are consistent with observations of the most powerful tidal disruption 
events, e.g.. Swift J1644. The jet velocity is, however, too low to match the Lorentz factor 
7 > 2 inferred in J1644. There is no such conflict in the case of other tidal disruption events. 
Since favorably oriented observers see isotropic equivalent luminosities that are highly super- 
Eddington, the simulated models can explain observations of ultra-luminous X-ray sources, at 
least in terms of luminosity and energetics, without requiring intermediate mass black holes. 
The spectrum remains to be worked out. Finally, since the simulated jets are bary on-loaded 
and have mildly relativistic velocities, they match well the jets observed in SS433. The latter 
are, however, more collimated than the simulated jets. This suggests that, even if magnetic 
fields are not important for acceleration, they may perhaps still play a role in confining the jet. 

Key words: accretion, accretion discs - black hole physics - relativistic processes - methods: 
numerical - galaxies: jets — X-rays: individual: Sw J1644-I-57 — X-rays: individual: SS433 
— X-rays: ULX 


1 INTRODUCTION 

Accretion flows around compact objects are commonly associated 
with relativistic jets. Such jets in active galactic nuclei. X-ray bina¬ 
ries in the low/hard state and gamma-ray bursts can be extremely 
powerful and can reach high Lorentz factors. It is believed that jets 
are driven by the rotational energy of a spinning black hole (BH) 
extracted via magnetic flelds plandford & Znajek|1977| BZ). The 
efficiency of energy extraction depends on the BH spin and the 
magnetic flux threading the BH horizon. The maximal power for 
a given spin is reached when the magnetic flux saturates at the 
value characteristic of the magnetically arrested disk (MAD) state 
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accretion flow. For rapidly rotating B Hs, the power of relativistic 
jets in the MAD state can exceed Mc^ ( [Tchekhovskoy et al.||2011 


jTchekhovskoy & McKinney|2012| t, i.e., more energy can go into 
the jet than the entire rest mass energy accreted by the BH. Such 
relativistic and magnetic jets have been extensively studied in re 
cent years, both numerically (e.g., | Tchekhovskoy et al.|20ffi |201 1 


Tchekhovsko y & McKinn ey |20 12 1 and analytically (e.g., |Penna et| 
al.|2013b[ Lasota et al. 12014] !. Although these studies have focused 


predominantly on hot, optically thin, geometrically thick accretion 
flows, recent numerical simulations ( |S^dowski et al.|2014a||McIG] 
[inney et al.|2014| | suggest that BH rotation-driven jets may be pro¬ 
duced with comparable efficiency even in the case of optically thick 
accretion flows. 


Relativistic jets are often invoked to explain the observed X- 
ray luminosities of tidal disruption events (TDEs). When a star is 
disrupted by a supermassive BH, part of the debris returns to the 
BH and is accreted at a very high, significantly super-Eddington 
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(“super-critical”) mass accretion rate. The isot ropic equivalent X- 
ray emission can reach as much as 10"^^ erg s"^ ( [Bloom et al.]2011 
[Burrows et al.|20lT| l. It is believed that this radiation is produced 
in internal shocks within a relativistic jet. However, as mentioned 
above, the most powerful jets require accumulation of a significant 
amount of magnetic fiux, and the disrupted star is incapable of pro¬ 
viding this fiux ( [Tchekhovskoy et al.|2014| l. In addition, a rapidly 
spinning BH is required to produce an effective Blandford-Znajek 
(BZ)jet. 

The initial mass accretion rate in a TDE is super-Eddington, 
which is very different from the highly sub-Eddington, radiatively 
inefficient regime that is usually considered when studying mag¬ 
netic jets. A super-Eddington accretion disk is dominated by radi¬ 
ation pressure and the gas is optically thick. It was predicted al¬ 
ready by [Paczynski & Wiita] ( [1980[ ) that such disks will be super- 
luminous. What is more, the large thickness of the disk naturally 
collimates radiation along the rotation axis, producing a signifi¬ 
cantly super-Eddington flux of energy in the funnel ( [Sikora|1981[ 
[Narayan et al.|198^ which can accel erate gas to mildly relativis¬ 
tic velocities ( Sikora & Wilson|1981| t. Such radiatively driven jets 
should occur in any super-critical accretion flow and should be 
present even if the BH is non-rotating or if there is little magnetic 
flux threading the BH horizon. In this work we study in detail the 
properties of such purely radiative (non-BZ) jets by means of gen¬ 
eral relativistic (GR) radiation MHD simulations. 

The paper is organized as follows. In Section[^we describe the 
numerical method. In Section[^we present the results, and discuss 
in detail the location of the photosphere (Section [TT] ), the proper¬ 
ties of the radiative jet ( [3.2[), the corresponding luminosities ( [3.3[ ), 
and the beaming factor In Section[^we discuss applications 
to TDEs, ultraluminous X-ray sources, and SS433. 


2 NUMERICAL METHODS 

We performed a set of three simulations of super-Eddington accre¬ 
tion on a non-rotating, a* = 0, moderately-supermassive, Mbh = 

3 X IO^Mq, black hole. We used the general relativisic (GR) ra¬ 
diation magnetohydrodynamical (RMHD) code KORAL ( |S^dowski[ 
[et al.|2013a[ [2014a[ ), which evolves gas and magnetic field in the 
ideal MHD approximation, together with the radiation field under 
the Ml closure scheme ( [Levermorej 19 84| ) . The models were initial¬ 
ized as optically thick equilibrium tori of gas and radiation in local 
thermal equilibrium, threaded by multiple loops of weak magnetic 
field. The only difference in the intitial states of the three simula¬ 
tions was the torus entropy parameter TC (see [Penna et al.|2013a[ l, 
which controls the initial torus density and thereby the mass accre¬ 
tion rate. 

The simulations were axisymmetric and were run with a reso¬ 
lution of 304 X 192 cells in radius and polar angle, with polar cells 
concentrated towards the equatorial plane. Each simulations was 
evolved for an exceptionally long time — roughly 2,000 orbital 
times at the innermost stable circular orbit (ISCO). This was pos¬ 
sible because the poloidal magnetic field was amplified and main¬ 
tained at a physically reasonable level by means of a mean-field 
dynamo, as described in [Sqdowski et al.[ ( [2015[ ). (Without the dy¬ 
namo, the poloidal magnetic field would quickly decay because of 
axisymmetry of these 2D simulations. On the other hand, without 
restricting ourselves to 2D, it would be virtually impossible to run 
the simulations for such long durations.) The duration of the sim¬ 
ulations was sufficiently long for the jet to reach equlibrium at the 


Table 1. Model parameters 


Name 

W 

tmax/ {GMIC^) 

<M>/MEdd 

1 

A 

10.0 

200,000 

45 

5.3% 

B 

5.0 

170,000 

310 

4.0% 

C 

1.0 

190,000 

4800 

4.6% 


Other parameters: BH mass: Mbh = 3 x 10^M©, BH spin: = 0.0, reso¬ 
lution: 304x192, radial extent of the domain: /?min = 1-85 to Rm&x = 10000, 
grid parameters (defined in [Sqdowski et al.|2015) Rq = 1.0, Hq = 0.6, 
minimal initial ratio of magnetic to total pressures: /5max = 10.0. 

outer boundary of the grid, located at r = 1000C[^ We used the 
following formulae for the absorption and scattering opacities, /Cabs 
and ATes, respectively, 

^abs = 6.4 X cm^/g, (1) 

/^es = 0.34cm^/g. (2) 

The model parameters of the three simulations are given in Table 
and are defined as in [S^dowski et al.[ ( [2015| ). 

The average accretion rates for models A, B and C were 45, 310 
and 4800MEdd. respectively, where the critical accretion rate Msdd 
is defined as, 

MEdd=^=2.44xlO‘«^gs-‘, (3) 

rjQC^ Mq 

where = 1.25 x lO^^M/M© erg s~^ is the Eddington luminosity, 
and 7/0 is the radiative efficiency of a thin disk for zero BH spin: 
rjo = 0.057. Eor the chosen BH mass (Mbh = 3 X lO^M©), L^dd = 
3.75 X 10^3 erg s"! and Mgdd = 7.32 x 10^^ g s-k 

Given the long durations of the runs (Table[^, the simulations 
reached inflow/outflow equilibriunj^ up to radius ~ 50 in the 
bulk of the disk, where the gas moves inward on average with a 
relatively low viscosity-driven velocity. In the region near the axis, 
where the gas flows out with quasi-relativistic velocities (as will 
be discussed below), the converged region extends out to the outer 
boundary at r = 10000. 


3 RESULTS 
3.1 Photosphere 

The simulated accretion flows are characterized by high, super¬ 
critical accretion rates and form geometrically thick disks (as 
expected). Eigure shows the time-averaged gas distribution in 
Model A, with color intensity indicating the mean density. Within 
the inflow/outflow equilibrium radius (req ~ 50) the disk den¬ 
sity scale height corresponds to hir ~ 0.3, a typical value for 
super-critical and radiatively inefficient accretion flows. Although 
most of the gas is concentrated near the equatorial plane, a sig¬ 
nificant amount is blown away to form the wind region visible 
in Eig. at intermediate polar angles. The lowest gas density is 
found near the polar axis. Eor Model A (accreting at the lowest rate, 
M - 45MEdd). the density on the axis at radius r = 1000 equals 
roughly 10"^^gcm“^. This corresponds to an optical depth against 
scattering greater than unity over the distance d = 1000 ^ 10^^cm. 
Thus, a significant fraction of the polar region is optically thick. 

^ Throughout the paper we adopt GM/c^ as the unit of length, and GM/c^ 
as the unit of time. 

^ Defined as a region where the local average radial velocity satisfies 
\v^\ > r/2t for a simulation of duration t. 
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Figure 1. Photosphere profiles for models A (black), B (red) and C (blue 
lines). Only model A has its photosphere inside the computational domain. 
The other two lines represent only the formal location of the photosphere 
which coincides with the domain boundary at r = 10000. The color levels 
represent the logarithm of density for model A. 


We estimate the radius of the photosphere in the funnel as a 
function of polar angle 6 by integrating the optical depth from the 
outer boundary of the simulation box down towards the BH along 
fixed 6. In computing the optical depth, we allow for relativistic 
effects and follow the formalism given in |S^dowski et al.| ( |2015| ) 
(their eqs. 47-49). The black contours in Fig. show the location 
of the photosphere for model A. Even for this moderately super¬ 
critical model, the photosphere at the axis is located as far out as 
r = 2000. For larger accretion rates the photosphere moves even 
farther out. In fact, for Models B and C, it lies otuside the compu¬ 
tational domain (red and blue lines in Fig. show only the formal 
photospheres coinciding with the outer edge of the domain). 


3.2 Radiative jet 

The left panel of Fig. shows again the density distribution in 
Model A, but limited now to the inner regions and with gas veloc¬ 
ity vectors superimposed. The right panel shows the magnitude and 
direction of the energy flux. To explain the quantity that is plotted, 
we start with the stress-energy tensors, Ty, Ry, of the magnetized 
gas and radiation, and note that -(T[ + R\) at any point is the radial 
flux of total energy (gas, magnetic and radiation) at that point as 
measured by an observer at infinity (lab frame). Also, pu\ where p 
is the gas density in the fluid frame and is the fluid four-velocity 
in the lab frame, is the radial flux of rest mass energy. Thus, 

r ^-{T^+R^^+pu^) (4) 

describes the flux of “useful energy” (i.e., rest mass energy sub¬ 
tracted out) as measured at infinity. Multiplying by then 
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Figure 2. Time-averaged results for model A in the poloidal plane. Left 
Panel: Color shows the distribution of the logarithm of density (in g cm“^) 
and arrows show the gas velocity. The white contour shows the formal bor¬ 
der of the jet region. Right Panel: Color shows the logarithm of isotro^ 
equivalent luminosity of total energy (in ergs“^ for definition see Eq.|^ 
and arrows show the direction of the energy flux. The white contour shows 
the border. 


gives the equivalent isotropic luminosity (ergs ^) that one infers 
from the measured flux at the given location, 

Liso - = -Anr^iT^ + + pw^). (5) 

The colors in the right panel of Fig. show this quantity, Liso, as 
a function of position, with arrows indicating the direction of the 
energy flux in the poloidal plan^ Over the entire volume, energy 
flows out, i.e., Liso > 0- In the region of the disk, bound gas is ac¬ 
creted, while in the wind and jet regions, unbound gas is ejected. 
Both correspond to positive Liso- The highest energy flux is in the 
region near the polar axis. Interestingly, even for the relatively 
low accretion rate of 45MEdd in model A, the isotropic luminosity 
reaches 5 x 10"^^ erg/s ^ ISOLEdd- 

To define the region corresponding to the jet, we follow 
[S^dowski et al.| ( |2013b| ) and require the local Bernoulli parameter 
p, defined as 

r^+Rr+pur ^ ( 6 ) 

pu^ pu^ 

to exceed 0.05, which corresponds to the velocity at infinity ex¬ 
ceeding roughly 30% of the speed of light. This definition is not 

^ The direction in the poloidal plane is determined by the vector (F^, F^) 
built from ortonormal radial and polar components of the energy flux as 
deflned in Eq.j^ 
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Figure 3. Top panel: Fractional contribution in model A of radiative {R\), 
magnetic (b^u^Uf - b^bf) and kinetic plus gravitational {pu^Ut + pu^) com¬ 
ponents of the energy flux to the total flux in the jet region as a function of 
distance. The shaded zone corresponds to the optically thin region above the 
photosphere. Bottom panel: Radial profiles of gas velocity in the jet region 
in model A. The solid line shows the gas velocity in the cell nearest to the 
axis, while the dotted line shows the density-weighted average velocity in 
the jet region. 


fundamental for our analysis, but is useful when talking about av¬ 
eraged or integrated quantities in the funnel. 

The white contour in Fig. [^denotes the border of the jet region 
defined according to the criterion given above. It agrees well with 
the region of lowest density (shown in the left panel) as well as the 
region of highest energy flux (right panel). 

The top panel of Fig. shows radial profiles of the fractional 
contributions of various components of the energy flux integrated 
over solid angle within the jet region. Red, purple and green lines 
correspond to the fluxes of radiative, kinetic, and magnetic ener¬ 
gies, respectively. The flux of thermal energy is negligible and is 
not shown. 

In the region of the jet closest to the BH (r < 1000), the energy 
flux is dominated by radiation which comes from the innermost and 
hottest part of the disk. However, the relative contribution of the 
radiative flux decreases with increasing distance from the BH as 
radiative energy is converted into kinetic energy (radiation pressure 
pushes the gas out). The magnetic contribution is unimportant at all 
radii. This is by construction because we specifically chose non- 
BZ-like conditions for these simulations - non-spinning BH, weak 
magnetic field well below the MAD limit. 

The transfer of energy from radiation to gas causes the gas 
velocity to accelerate along the funnel, as seen in the bottom panel 
of Fig. The gas at the axis reaches roughly 50% of the speed of 
light at distance r ~ 1000. The density-weighted gas velocity in the 
funnel reaches only ~ 35% of the speed of light, reflecting the fact 
that gas near the edge of the funnel is denser and less effectively 
accelerated. 

Acceleration continues as long as momentum is transferred 
from the radiation, i.e., up to radius r ^ 2000. Constant gas ve¬ 
locity implies that the kinetic energy has saturated. However, the 



Figure 4. Radial velocity of the gas at the axis (solid lines) and the maximal 
velocity the gas will achieve if all the energy is converted into kinetic energy 
at infinity (dashed lines). Results are shown for all three models. 

fractional contribution of the kinetic energy (purple line in the top 
panel) increases even further. This is because of the residual radi¬ 
ation crossing the borders of the jet region and increasing the ki¬ 
netic contribution. Acceleration stops when radiation drag by pho¬ 
tons from the funnel walls is large enough to counter-balance the 
acceleration by radiation flowing along the funnel ( [Sikora & Wil-| 
|son In the framework of the Ml closure scheme, this cor¬ 

responds to the situation where the gas comoves with the radiation 
rest frame, where no momentum is transferred between the radia¬ 
tion and the gas. 

Even though we ran our simulations over a large domain ex¬ 
tending to Tout = 10000, it turned out that the outer boundary was 
still not far enough out to resolve the photospheres of Models B 
and C. In these models, gas acceleration occurs all the way to the 
boundary. Nevertheless, we can still estimate the terminal velocity 
of the gas by looking at the total energy flux near the boundary 
and assuming that all the energy is converted into kinetic energy 
at infinity. The velocities along the axis and the terminal velocities 
for the three models are shown in Fig.j^by solid and dashed lines, 
respectively. The gas in the lowest accretion rate model A reaches 
saturated velocity 0.5c already at radius r - 1000. The terminal 
velocity for this model (black dashed line) is only ~ 10% higher, 
which means that almost all the available energy has been converted 
into kinetic energy. For the other two models, with accretion rates 
310 and 4800MEdd. the actual and terminal gas velocities are signif¬ 
icantly lower and reach only 20 to 30% of the speed of light. This 
is probably the result of more effective mass loading of the funnel 
region in these models. Because of the larger optical depth, we pre¬ 
sume radiation more effectively drives gas towards the axis, which 
results in weaker acceleration of the gas in this region. 

3.3 Luminosities 

Fig.|^presents the equivalent isotropic luminosity (top panel, Eq.[^ 
and the radial velocity of the gas (bottom panel) as a function 
of polar angle for models A (M = 45MEdd). B (310MEdd), and C 
(4800MEdd)- The luminosities were measured at radius r = 5000 
and reflect the flux of total energy, which in the funnel/jet region 
is dominated by kinetic luminosity. The solid lines correspond to 
the region of the simulation domain where inflow/outflow equilib¬ 
rium state has been reached at r = 5000. The dotted lines refer to 
unconverged region, which is not in flow contact with the BH and 
therefore the results are less reliable. 
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Figure 5. Top panel: Isotropic equivalent luminosity (measured at r = 
5000) for models A (black), B (red) and C (blue lines) as a function of the 
polar angle 6. The dashed lines show the fitting function given in eq. □for 
a black hole mass M = 3 x lO^M©. Bottom panel: Velocity of the gas at 
r = 5000. In both panels the solid lines correspond to the region of the sim¬ 
ulation that is in infiow/outfiow equilibrium at r = 5000, while the dotted 
lines refer to unconverged regions. 


For an observer looking precisely down the axis, Liso = 5 x 
10"^^erg s"^ for model A, and increases to as much as 3 x lO'^^erg s~^ 
for model C (4800MEdd)- The isotropic equivalent luminosity is 
lower for off-axis observers, dropping by roughly an order of mag¬ 
nitude between the axis and the border of the jet region. The lu¬ 
minosities scale linearly with accretion rate, i.e., the efficiency of 
energy injection into the funnel region seems to be independent of 
the accretion rate. 

The isotropic equivalent luminosity profiles are reasonably 
well described by the following fitting function. 


Lfl, = 4 X 10'^^ 


lOWfidd 


(7) 


where the scaling with BH mass is includecQand where ^jet,0 - 0-2 
is the characteristic jet beaming angl^ 

The bottom panel of Fig. shows the dependence of the gas 
velocity on polar angle for each of the models. Model A with the 
lowest accretion rate shows the most efficient acceleration. In all 
three models, the velocity quickly decreases with increasing 6 and 
transitions to < 0.1 in the wind region (6 > 0.5). 


^ Although we simulated accretion flows on a supermassive BH (to directly 
show how strong the radiative jet can be) our results can be applied to any 
astrophysical BH by scaling them properly with mass. This is because the 
optical depth over given number of gravitational radii, which determines, 
e.g., the efficiency of gas/radiation interaction, is independent of BH mass 
for fixed fraction of Eddington accretion rate. 

^ Allowing for the fact that there is more solid angle at larger 6, a better 
estimate of the luminosity beaming angle is twice this: 0jet ~ 0.4. 
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3.4 Beaming factor 


Disks accreting at super-critical rates are geometrically very thick. 
As a result, energy in general and photons in particular escape pref¬ 
erentially along the axis, where the density and optical depth are 
lowest. Most of the energy escapes through the funnel, and the 
flux of energy decreases steeply with increasing polar angle (Eq.[7] 
Fig.0. It is useful to define the beaming factor b(6), which gives 
the ratio of the isotropic equivalent luminosity Liso (defined in eq.[^ 
inferred by an observer at polar angle 6 to the real luminosity of the 
source Ltot, 

J Tiso Tiso 

Ltot TjMc^ 

where Ltot is obtained by integrating the energy flux over the 
spher^ 

Ltot = dM<p = rjMd, (9) 

with rj being the total efficiency of the accretion flow (results given 
in Table [^. When integrated over the sphere, the beaming factor 
should, on average, be equal to unity. However, because the equi¬ 
librium region at r = 5000 is limited to the polar region (gas near 
the equatorial plane was not able to relax yet because of long vis¬ 
cous timescale there), this integral cannot be directly calculated. 

Fig. shows the beaming factor as a function of 6, as mea¬ 
sured at radius r = 5000. For all three runs, b is maximum at the 
axis and decreases away from the axis. The maximal beaming fac¬ 
tor is Z? = 3.5 for model A, and b = 2.0 - 2.5 for models B and 
C. More pronounced beaming in model A may come from the fact 
that, in this case, b is measured outside the photosphere. The beam¬ 
ing factor equals unity, i.e., local isotropic equivalent luminosity is 
equal to the true luminosity, at ^ ^ 0.35, which roughly coincides 
with the edge of the funnel. In the wind and disk regions, b{6) 1, 

reflecting the fact that comparable amounts of energy leave in the 
funnel and in these regions, but the solid angle covered by the fun¬ 
nel is a small fraction of the whole sphere. 

It should be emphasized that we defined the beaming fac¬ 
tor b(6) with respect to the total luminosity of the accretion flow, 
not the Eddington luminosity. The three simulations discussed in 
this paper are all energetically efficient, as reflected by the fact 
that they have efficiencies 77 (Table nearly equal to the standard 
Novikov-Thorne thin disk efficiency, 770 = 0.057. Correspondingly, 
Ttot/TEdd ~ MI^Edd- Therefore, if we define an Eddington-scaled 
beaming factor we roughly expect 


^Edd ^ 


Tiso 

^Edd 


J^b. 

^Edd 


( 10 ) 


Considering that b itself is of order a few along the axis, it is clear 
that Z?Edd can be extremely large for highly super-Eddington accre¬ 
tion flows. If the BH was spinning and if it was threaded by a MAD- 
level magnetic field, the luminosity would be larger still by possibly 
a factor of several. 


4 DISCUSSION 

Magnetically-dominated, BH rotation-powered jets in radiatively 
inefficient hot accretion flows have been simulated and studied now 

^ In steady state, the integrated energy flux is independent of radius, so it 
can be evaluated at any convenient radius. The values of Ltot and 77 quoted 
here correspond to r = 5000. 
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Figure 6. Beaming factor (for definition see Eq. as a function of polar 
angle in the three simulations. 



for a number of years. These jets operate through the BZ mecha¬ 
nism, in which ordered magnetic fields are twisted by the rotating 
BH. The twisted field lines unwind and expand in a direction par¬ 
allel to the BH spin axis because of the collimating action of the 
geometrically thick accretion disk. 

The Jets studied in this paper involve different physics. Our 
BHs do not rotate and there is no strong ordered magnetic field, so 
there is no opportunity for the BZ mechanism to operate. Neverthe¬ 
less, the super-Eddington accretion flows considered here do pro¬ 
duce powerful jets. Because these accretion disks are geometrically 
thick, they naturally form funnels oriented parallel to the angular 
momentum axis of the accretion flow (just as in the case of radia- 
tively inefficient hot accretion flows). MRI-driven turbulence in the 
disk leads to energy dissipation in the disk interior, which heats the 
gas. The gas, in turn, cools by emitting radiation. The radiation then 
flows down the gradient of density/optical depth via diffusion and 
turbulent transport ( [Jiang et al.|2014| |, which naturally concentrates 
a significant amount of energy in the funnel region. There, radia¬ 
tion pressure accelerates gas outward along the axis of the funnel, 
resulting in a jet that travels at a modest fraction (0.2-0.5) of the 
speed of light. 

The process described above will take place so long as the 
accretion flow forms a thick disk with a funnel. Independently of 
whether the funnel is optically thin or thick, radiation always moves 
towards the axis. Once there, it collimates and rushes out along the 
axis, transfering its momentum to the gas. The transfer operates as 
long as there is enough optical depth. Once the radiation crosses the 
photosphere, where radiative and kinetic energies are of the same 
order (see Fig.[^, the acceleration of the gas terminates. 

Efficient conversion of radiative to kinetic energy happens 
only if there is sufficient optical depth along the funnel. Already 
for an accretion rate of 45MEdd (model A, the lowest M we have 
considered), most of the radiative energy is converted to kinetic en¬ 
ergy. Eor lower accretion rates, the optical depth in the funnel will 
be lower, and we expect that a larger fraction of the jet energy will 
remain in the form of radiation rather than gas kinetic energy. How¬ 
ever, the total jet power should scale similarly for all super-critical 
accretion rates. This is because the efficiency of accretion and the 
thickness of the disk (which collimates the jet) are independent of 
the accretion rate in this regime ( jS^dowski et al.|2015] ). 

Can the kinetic energy escaping the funnel be observed by a 
distant observer? Yes, provided it is converted back into radiation. 
In magnetic jets, reconnection is possibly a dominant mechanism 
for this conversion ( jSironi et al.|20I5] ). However, magnetic recon¬ 
nection is unlikely to be important for the radiatively-driven jets 
considered here since there is so little magnetic energy (Eig. 3). 
Shocks are a more natural mechanism for converting the gas ki¬ 


Figure 7. Velocity profiles along the axis for model A corresponding to three 
epochs. 


netic energy to radiation via nonthermal particle acceleration. To 
show that shocks are quite likely. Fig. [^presents the radial profile 
of gas velocity near the axis for model A at three different moments 
of time. We see that the velocity profiles along the funnel are highly 
variable. Layers (shells) of gas moving with higher than average 
velocity will clearly overtake preceding lower-velocity gas, thereby 
producing shocks (especially in optically thin regions). Such “inter¬ 
nal shocks” will dissipate a good fraction of the kinetic energy, con¬ 
verting it to thermal and nonthermal energy and ultimately into ra¬ 
diation. What fraction of the total energy flux measured at r = 5000 
(Fig.[^ reaches the observer in the form of radiation depends on the 
geometry and dynamics of this process. The energy distribution of 
the reprocessed radiation can be significantly different from that of 
the original radiation in the funnel. Estimating the spectral shape 
is, however, beyond the scope of this paper. 


4.1 Comparison with previous works 

Studies of geometrically thick disks corresponding to high accre¬ 
tion rates were initiated in the late 1970s by authors studying thick 
equili brium torii ( |Fish^ne^^^lo ncrief| 1976[ [Kozlowski et al.j 


|1978||Abramowicz et al.|1978||Jaroszyns]d et al.|1980j l. It was re- 

alised that if torii are thick and radiative pressure supported, their 
total luminosity must exceed the Eddington limit ( jPaczynski & Wi-j 
|ita|1980j). Later w orks by |Sikora| ( |1981j ), |Sikora & Wilson| ( |1981j ) 
and jNarayan et al.| ( |1983j ) specifically analysed the properties of the 
polar outflows. The accretion flows were modeled by thick equilib¬ 
rium torii with arbitrarily prescribed angular momentum profiles. 
It was shown that, if such torii are thick enough, the local radiative 
flux in the funnel may significantly exceed the Eddington value. 
jSikora & Wilsonj ( |1981| ) showed in addition, that the flow cannot 
convert more than a small fraction of the disk’s radiative luminosity 
into kinetic energy in the polar outflow because of radiation drag. 
However, the authors noted that the acceleration may be more effi¬ 
cient if the funnel is optically thick. These predictions are in perfect 
agreement with our study. 

Numerical simulations of super-critical accretion are rela¬ 
tively new. The field was initiated in a series of seminal works by 
jOhsuga et al.j ( [2005 [ [2009 j ); jOhsuga & Mineshi^ ( [2007 j [201 Ij i, 
who performed radiation hydrodynamical, and later MHD, simula¬ 
tions using a non-relativistic code and adopting the flux-limited dif¬ 
fusion approximation. More recently, our group developed general 
relativistic numerical codes with Ml radiative closure and applied it 
to super-critical accretion ( jS^dowski et al.|2014a||2015[|McKinney| 
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let al.|2014> . A similar scheme was implemented in an independent 
radiation GRMHD code by [Fragile et al.| ( |2014j ). Most recently, 
[Jiang et ^ ( |2014j ) simulated a super-critical accretion flow using 
a Newtonian code with cylindrical coordinates (and a cylindrical 
BH). The authors solved the radiative transfer equation directly on 
a flxed grid of angles, instead of applying any radiative closure. 

All the studies cited above agree qualitatively on the main 
features of super-critical accretion, viz., that there is no barrier to 
super-Eddington accretion, that these flows are optically and ge¬ 
ometrically thick, and that they are characterized by strong polar 
outflows and signiflcant winds. The detailed composition of the 
outflows differs between studies, e.g., in [Jiang et al.[ ( [2014] ) energy 
flows out primarily via radiation, whereas our jets are kinetic en¬ 
ergy dominated, and the reason for this is presently unclear. Our 
work is the first one to study polar outflows using a general rela¬ 
tivistic code. Also, we simulated accretion onto a (mildly) super- 
massive BH (all previous work considered stellar-mass BHs). 

In a work similar to ours, [Takeuchi et alT| ( [2010[ ) studied radia- 
tively driven jets using a Newtonian code and simulated an accre¬ 
tion flow with M ^ 6MEdd- The authors focused on the acceleration 
and collimation mechanisms of the radiative jet. They showed that 
the jet can be accelerated up to mildly relativistic speeds by the 
radiation-pressure force and collimated by the Lorentz force of a 
magnetic tower structure in the innermost region. They also pointed 
out the importantce of the radiative drag force in limiting gas ac¬ 
celeration in [Takahashi & Ohsuga[ ( [2015] ). The gas velocities they 
obtained are in perfect agreement with our study. However, we do 
not observe magnetic collimation of the outflow. This discrepancy 
may be related to the fact that the simulated flows in the present 
work have significantly higher accretion rates. 

4.2 Astrophysical implications 

A principal result of this paper is that super-critical accretion on 
BHs leads to a collimated outflow of energy along the polar axis 
even if the BH spin is zero and BZ-like extraction of BH rotational 
energy is ruled out. For observers looking directly down the funnel, 
the isotropic-equivalent luminosity can significantly exceed the Ed¬ 
dington luminosity (Eq.|7]). The observed luminosity drops with in¬ 
creasing polar angle. We discuss the importance of these results in 
the context of tidal disruption events (TDEs), ultra-luminous X-ray 
sources (ULXs), and the microquasar SS433. 


suggestion that the extreme luminosity was produced by a relativis¬ 
tic jet driven by rotational energy of the BH. However, efficient ex¬ 
traction of BH spin energy requires the accumulation of signiflcant 
magnetic flux at the BH, far more than the disrupted star can pro¬ 
vide ( [Tchekhovskoy et al.|2014[ |. This problem could be overcome 
if the tidal stream is able to drag magnetic held from a pre-existing 
accretion flow around the BH ( [Kelley et al.[2014[ ). 

The initial highly super-Eddington regime of accretion in a 
TDE corresponds precisely to the simulations discussed in this 
study, where we consider super-critical accretion on a supermas- 
sive BH with M of order tens to thousands of Eddington. We can 
therefore explore whether a purely radiative jet, similar to the ones 
studied in this work, can explain observations of TDE jets. As 
shown in Section [T^ for an observer looking down the funnel, 
the isotropic equivalent luminosity reaches (Eq. 0 10"^^ergs ^ for 
canonical TDE parameters. By increasing the BH mass by an order 
of magnitude to lO^M©, it is possible to match the peak luminos¬ 
ity observed in Swift J1644. However, the radio emission in the 
afterglow of this TDE event suggests that the jet must have been 
relativistic, with Lorentz factor > 2 ( [Berger et al.[2012[ |. Such high 
velocities are difficult with purely radiative acceleration — the ra¬ 
diation held in the funnel is not sufficiently collimated, and acceler¬ 
ation is prevented by radiative drag. Therefore, it seems that radia¬ 
tive jets alone cannot explain the properties of Swift J1644. Instead, 
the jet in this object probably involved a combination of radiatively 
and BH rotational energy driven acceleration. 

Some other TDE candidate events, e.g., PSl-lOjh ([Gezari et 
al.|2012>, P Sl-llaf ( |Chomock et al.|2014> , and PTF09ge lAr^ 
et al.||2014) , did not show significant X-ray and radio emission, 
but were observed in the ultraviolet band only. This fact may be 
explained if we assume that these sources do not produce strong 
magnetic jets accelerating gas to such high velocities as in the 
case of Swift J1644. Having in mind that the onset of accretion 
of the returning debris corresponds to strongly super-critical accre¬ 
tion, and that the luminosities are super-Eddington, the observed 
optical emission would come from the funnel region where it is 
reprocessed through the kinetic flux, as discussed in this work. It 
could also come directly from the inner funnel wall provided the 
optical depth in the polar region is low, though this is likely only 
at lower accretion rates. In both cases, it is required that the TDE 
candidates must be observed within the jet beaming angle, because 
the luminosity would be significantly lower for more inclined ob¬ 
servers (Fig. 


4.2.1 Tidal disruption events 


4.2.2 Ultraluminous X-ray sources 


Whenever a star comes closer to a supermassive BH than its tidal 
radius, it is disrupted by the tidal gravitational held ( [Rees[p^88l 
[Evans & Kochanek|1989[ t. Approximately half the mass of the star 
remains bound to the BH and returns after a while to the vicinity of 
the BH. Relativistic precession causes the tidal stream to interact 
with itself, shock and circularize (e.g., [Bonnerot et al.|20T5l[Sh-[ 
[iokawa et al.|2015| l. Once the gas circularizes, it starts to build up 
magnetic held which triggers accretion. According to the standard 
model (e.g., [Krolik & Piran|2012[ ), for typical parameters (stellar 
mass star disrupted by a ~ lO^M© BH), the peak accretion rate at 
the onset of accretion can exceed 10^M^dd- This accretion rate sub¬ 
sequently declines as 

In the case of the TDE Swift JI644, the observed isotropic- 
equivalent X-ray luminosity was as high as lO'^^erg s"^ during the 
initial outburst peak ( [Bloom et al.[2011[ [Burrows et al.[201l| l. This 
fact, together with the properties of the radio-emission, led to the 


Ultraluminous X-ray sources (ULXs) are objects emitting isotropic 


equivalent luminosities in excess of 10^^ erg/s (e.g.. Miller et al. 


[2004[ ). The critical number corresponds roughly to the Eddington 
luminosity for a lOM© central source. A standard thin disk around 
a IOMq BH, accreting at a sub-Eddington rate, cannot produce the 
luminosity observed in ULXs, so one possibility is that the accret¬ 
ing compact objects are intermediate-mass BHs. However, if the 
emitted radiation is beamed into a relatively narrow solid angle, the 
observed luminosity can easily exceed the Eddington value for ob¬ 
servers located in favorable directions, even though the total emis¬ 
sion may still be sub-critical. In this scenario, ULXs are normal 
lOMo BHs ( [King|2009[ ). 

Our study supports the latter picture since we And very strong 
beaming of radiation in geometrically thick disks. The particular 
systems we studied have very high and extreme accretion rates. In 
these flows the funnel region is optically thick and radiative energy 
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is efficiently converted into kinetic energy of the outflowing, col¬ 
limated gas. Whether or not this kinetic energy can be dissipated 
in shocks to produce radiation and whether the resulting spectrum 
would match observations of ULXs is uncertain. 

However, another regime of super-critical accretion may be 
more applicable to ULXs. For not-so-extreme accretion rates, say 
M ~ lOMEdd. the polar region is optically thin and conversion of 
radiative to kinetic energy is not efficient. It is then possible to ob¬ 
serve directly the radiation emitted inside the funnel. Although the 
beaming may be less strong, it can still magnify the observed lumi¬ 
nosity and cause a stellar-mass BH source accreting at a moderately 
super-critical rate to match the isotropic-equivalent luminosities 
observed in ULXs. Detailed modeling of the observational prop¬ 
erties of such systems was done by |Kawashima et aL| ( |2012| ) using 
simulation data corresponding to a mildly super-critical (~ lOMEdd) 
accretion flow on a lOM© BH. The simulation was done with a 
Newtonian radiation MHD code. The authors showed that Comp- 
tonization in the polar region significantly hardens the original pho¬ 
ton spectrum. 

We have estimated the beaming factor, defined as the ratio of 
the observed (isotropic equivalent) to the real (integrated over the 
sphere) luminosities, and find (Fig.[^ that it significantly decreases 
with increasing polar angle. Therefore, one should expect that, for 
the same BH mass and accretion rate, the distribution of ULX lu¬ 
minosities will be sensitive to the orientation on the sky. 


4.2.3 SS433 


SS433 is the first microquasar to be discovered. It is an eclips¬ 
ing X-ray binary system in which the compact primary object is 
most likely a black hole with Mbh ~ lOM© ( [Grindlay et al.|1984] ). 
The wind of the companion star feeds the central object at an 
extremely high rate. The measured disk wi nd outflow rates are 
M = 10-5 _ JO--* Moyr‘ ^Shkl ovskylp^^ which suggests that 
the accretion rate on to the compact object is super-critical. A di¬ 
rect measurement of the accretion M is, however, impossible be¬ 
cause the innermost regions of the disk are obscured. SS433 shows 
powerful jets precessing with a period of 162 days and sweeping 
twin cones with half-angle 20°. The jets are loaded with baryons 
and move at a velocity of 0.26c ( |Margon|1984[|Fabrika|2004| l. The 
corresponding kinetic powe r dominates the jet en ergy budget and 
equals 10^^^ - 10^^-^ erg s"^ ( [Marshall et al.|2015 l. 

The fact that the jets are precessing suggests that the BH spin 
is non-zero and that there is a misalignment between the angular 
momentum of the binary and the BH spin axis. These conditions 
have to be satisfied for Lense-Thirring precession ( [Bardeen & Pet-j 
[terson|1975] ) to occur. If the BH is spinning, then it seems reason¬ 
able to assume that the observed jets are powered by the BH ro¬ 
tational energy via a strong magnetic field. However, the origin of 
the baryons in the jet is then problem since the magnetic field lines 
will screen the baryons in the disk from the jet region. It has been 
suggested that the jets are initially leptonic, but they then acquire 
baryons when they collide either with the thick accretion disk itself 
or the disk wind. 

We propose another explanation for the baryon-loaded jets in 
SS433. The mass supply rate on the compact object most likely cor¬ 
responds to a super-critical accretion. At such high accretion rates, 
an optically thick and geometrically thick disk will form and, as 
shown in this work, radiation moves into the polar region, push¬ 
ing gas along with it and making the funnel optically thick. The 
radiation then escapes up the funnel and accelerates the gas, thus 
converting a good fraction of the total energy flux into jet kinetic 


energy. Once such a radiatively-driven jet crosses the funnel photo¬ 
sphere, it is baryon-loaded, kinetically-dominated and collimated. 

These features in our simulations are consistent with observa¬ 
tions of the jets in SS433. The density-averaged gas velocity in our 
simulated jets is of the order of 0.2 - 0.3c (bottom panel of Fig. [^, 
which is in excellent agreement with the velocity observed in 
SS433. For a lOM© BH accreting at lOOMEdd, the jet power we ex¬ 
pect based on the simulations is roughly O.OIMc^ ^ 2x 10"^*^ erg s"^. 
This is more than adequate to explain the observed kinetic luminos¬ 
ity of the jets in SS433. 

The non-zero BH spin and misalignment between the orbital 
and BH spin axes, both of which are required for jet precessions, 
suggest that additional effects over and above simple radiative ac¬ 
celeration might be operating. The small observed opening angle of 
the jet in SS433 (1° - 2°), is not consistent with the jets seen in the 
simulations, which have opening angles of the order of 10° - 15° 
(Fig.§. Perhaps there are magnetic effects on top of the radiative 
processes we have discussed, which cause strong collimation. 
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